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Abstract. In order to develop a toy model for car's traffic in cities, in this paper we analyze, 
by means of numerical simulations, the transition among fluid regimes and a congested jammed 
phase of the flow of kinetically constrained hard spheres in planar random networks similar to 
urban roads . 

In order to explore as timescales as possible, at a microscopic level we implement an event driven 
dynamics as the infinite time limit of a class of already existing model {Follow the Leader) on 
an Erdos-Renyi two dimensional graph, the crossroads being accounted by standard Kirchoff 
density conservations. We define a dynamical order parameter as the ratio among the moving 
spheres versus the total number and by varying two control parameters (density of the spheres 
and coordination number of the network) we study the phase transition. 

At a mesoscopic level it respects an, again suitable adapted, version of the Lighthill-Whitham 
model, which belongs to the fluid-dynamical approach to the problem. 

At a macroscopic level the model seems to display a continuous transition from a fluid phase 
to a jammed phase when varying the density of the spheres (the amount of cars in a city-like 
scenario) and a discontinuous jump when varying the connectivity of the underlying network. 
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1 Introduction 



In the past decades an always increasing interest has been payed to the flow of cars in 
urban roads (see e.g. [1] or [2] for a beautiful modern review): a primary challenge is the 
reduction of time delays and CO emissions, in a nutshell, the congested traffic flow [3]. 
Despite progresses developed by engineers and physicist (see i.e. [HE]), essentially fo- 
cused on large streets or small tree-like graphs [6l [7] , very little is known concerning the 
behavior of traffic on large two dimensional networks [2]. 

In a completely different context, the last twenty years saw the statistical mechanics of 
disordered and complex systems [S] experience an increasing development as well as its 
range of applicability (see for instance [9], [TOl [11]) and, inspired by these successes, we 
want to investigate the nature of the transition among a fluid state and a jammed coun- 
terpart in traffic flow on planar networks, investigating the presence (or the lacking) of 
criticality fl2[ [T3] within a out of equilibrium statistical mechanics framework [T^ [T5] . 
Furthermore, statistical mechanics of disordered systems recently pointed out a deep con- 
nection among replica symmetry breaking scenario [T6l [8] , the paradigm of the transition 
among fluid and glass and the P —>■ NP transition in problem solving of hard satisflabil- 
ity problems [171 [H] • 

Interesting, if the mapping among jamming transition and P NP completeness would 
apply to traffic jams too, it would vanish every attempt to an online control of car flow 
by external massive macro-computing giving more flrm ground to interacting local opti- 
mizers (as i.e. neural networks) [T9l [20| [21]. 

Deepening our knowledge concerning the jam transition in traffic flow should be then of 
great importance, in our traffic optimization planning, if, varying tunable parameters, 
glassy-like criticality arises [22] . 

From a practical viewpoint, as a rigorous formulation of out of equilibrium statistical 
mechanics is far from being exhausted [231 [2H [H] , we do not have a paved mathematical 
way to follow for checking i.e. the involved time-scales [25] or the reach of a stationary 
state [2n] (which is a primary requirement for giving meaning to the averages) and con- 
sequently there is the need of fastest simulation algorithms [271 [2S] to cover as timescales 
as possible. 

Even though we will move toward a molecular-dynamics-like approach [29], we stress 
that for similar reasons the biggest amount of works on this subject uses in fact cellular 
automata [301 [SB [3 [32] which are quite faster than the continuous models [331 [Ml [35] . 
Fast simulations are in fact hard tasks, especially in models with continuum potentials 
as they need to be made discrete generally by using Trotter expansions of the Liouvillean 
[29] describing the motion in the phase space, forbidding a very long simulation time (i.e. 
by Liapounov constrictions [36]). 

Avoiding these potentials, a very fast integration of the dynamics is offered by the Verlet 
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event driven dynamics of hard spheres [37]: these spheres are without a real potential; 
they move on straight lines, up to a core distance at which they touch one-other and they 
feel an infinite barrier of potential which converts instantaneously the kinetic energy into 
potential energy. 

In this way the motion against two successive collisions does not require integration. It is 
in fact propagated from collision to collision, the new positions and momenta are worked 
out by imposing conservation of particles (Kirchoff rule), energy and momenta, (we will 
preserve just the Kirchoff rule in our framework) and the motion is propagated again and 
so on [3H] (we emphasize that this approach has been tackled also to granular systems 
[39] , which are glassy systems sharing several features with traffic flow [iQl SI] ) . 
Of course there exist already several very sharp models for traffic flow but we introduce 
our one because we are moving in an opposite way for a different scope with respect the 
standard approach: as we want a large amount of cars as well as long simulation time 
for a thermodynamical approach, we allow ourself to skip as details of the motion as 
possible, retaining just the main features (as usually happens when looking at criticality 
in statistical mechanics |42j). 

Once defined the microscopic dynamics we then focus at first at the mesoscopic level 
(order 10^ cars) to recover the Lighthill-Whitham scenario [13], for showing consistence 
with pre-existing works, then we focus on a macroscopic scale (up to 10"^ cars) to study its 
thermodynamics. We introduce the ratio of the moving particle as a standard dynamical 
order parameter [2], labeled by that we call fluidity for the sake of clearness, and define 
it as 

N 

(p=l-N~'J2v„ Vie[0,l] (1.1) 

i 

(such that it is trivially one in the jammed phase (where there is no longer any motion) 
and decreases toward zero in the liquid phases) and study its behavior: average, distri- 
bution and fluctuations. 

The model seems to display a jammed phase where the fluidity is strictly zero and its 
fluctuations are delta-like centered on the average (the congested phase where all the 
spheres are caged among their nearest neighbors) and a flowing phase in which the flu- 
idity seems to decrease continuously to zero (cages smoothly disappear) by decreasing 
the density or discontinuously by increasing the connectivity and its fluctuations appear 
Gaussians. 

The whole suggesting the model undergoes a second order like transition in the density 
and a first order like in the connectivity. 

For the sake of clearness we aim to label with a the connectivity of the network, with 
p the density of the cars (even thought, from practical comparison of different size 
networks it will be sometimes easier to deal directly with the un-normalized amount of 
car A^) and with Vi = v Wi E {1, N) the velocity of the i^'^ car. 
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It is worth noting that the standard technique of statistical mechanics on diluted systems 
[1UI15] merges the two parameters via the relation atanh(p) = p', p' being an equivalent 
density in a fully connected network, so actually do not display clearly the transition 
split among the two control parameters 

Furthermore we want to stress that there exist already works on dilute hard spheres in 
different contexts (as on the Bethe lattice) [IZlllH] but, to our knowledge, not on networks 
with topology close to urban one. 



2 Microscopic model 

In this section we point out the simulation scheme, which, for the sake of simplicity we 
spit in two parts: the choice of the underlying network (the topology) and the choice of 
the dynamics on the network (the interactions). 



2.1 The Network 

At first we must introduce the graph. In order to mimic a real urban center we think 
at a graph whose links represent the roads and vertices represent intersections and end 
points. For their high connectivity, scale free |19] and fully connected [50] networks are 
inappropriate to describe such a graph and for the extremum order and homogeneity 
they present, also Voronoi tessellation [51] and regular grids are avoided [52] . 
Real data show a linear dependence among the number of roads versus the number of 
intersections [53], whose ratio must be obviously between one (tree like structures with 
no loops) and two {2D regular lattice) with an empirical slope close to 3/2 [SI]. Even 
though clever growth algorithms for these networks recently developed [SI], for compu- 
tational simplicity (as we will have to average over several configurations and we want 
the fastest procedure) we choose the planar Erdos-Renyi graph [39] above the giant com- 
ponent threshold, which is close to the requested class of random graphs [55] and is of 
immediate realization on a computer as no growing is concerned. 

All the links represent streets built by two lanes so to have both an incoming and an 
outgoing flux from each node. 

On this network we can vary its averaged connectivity -the coordination number- (de- 
noted by a) so to explore from the region of extreme dilution near the percolation thresh- 
old up to a fully ordered grid. 

Of course to check convergence to the infinite volume limit we will test our simulations 
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Figure 2.1: Examples of generated networks. At left an high connectivity topology, at 
right a low one. We start from a fully ordered grid in which each node has exactly 8 
links then we dilute the links with a Poisson distribution with mean 4, such that the main 
nodes are ordinary crossroads. If a node has two, one or zero links is removed too. 



varying also the size of the grid. 



2.2 The dynamics 

On this network we let live N cars for which we want to menage two limits at the same 
time (which conflict in terms of CPU time): the large N limit to have enough data for 
the averages and the infinite time behavior so to approach to a steady state: the need of 
the fastest plausible algorithm follows straightforwardly. 

One of the most important class of models which aim to mimic car dynamics is the so 
called follow the leader [56l[57l[6], where the dynamics of the i*'' car is assumed to respect 
the following differential equation (or some suitable variants [2]), where {i + 1)*'' labels 
the car in front of the i*'', following the direction of motion: 

^-(^.--m). (2.1) 

In a nutshell the car i accelerates if the car in front of is accelerating too (this happens 
both for positive and negative accelerations). 

As a solution for the long time steady state behavior of the follow the leader model is 
given by the same constant velocity for all the cars (as it appears clearly by solving 
eq. fl2.ip in the large time limit by imposing d'^Xi/dt^ = 0) we choose for our microscopic 
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dynamics an event driven hard-sphere-like dynamics [38], which is known to be one of 
the fastest achievable dynamics in terms of CPU time [281 EH] and is in agreement with 
these continuous potential in the regime where we are interested in [t —>■ oo). 
On every lane overtaking is forbidden (FIFO principle [5Z]) and every car has an energy 
Ei = vf/2 during the free motion, whose dynamics is propagated between collisions 
among two cars or one car and a crossroad. 

There the collision rules are worked out with a remarkable difference with respect to 
canonical physics: nor detailed balance neither the third law of dynamics do hold. When 
a collision happens there is no conservation of momenta and the incoming car loses all 
the energy then starts again following the collided car (crossroads apart which confer 
a certain degree of randomness), for this reason we call our hard spheres kinetically 
constrained [59], [60]. However as for instance elegantly explained in [611 |T5] this is not 
a too serious limitation when investigating out of equilibrium steady state (the same is 
not true of course in equilibrium [12] ) . 

The potential felt by the cars is the classical hard core potential [321 [2B] V, which, by 
using \xij\ to evaluate the distance between the two cars situated respectively at Xi and 
X j , can be written as 

V{\x.,\)^l J' = ^ h|-^j (2.2) 

^' -"^ \ Vi = oo if |xij| < d j ^ ^ 

where d can be thought of as a security distance, the minimal distance allowed among 
two consecutive cars (we stress that varying d changes also the maximum number of 
admitted cars inside a network Nmax -which sets p = 1- and consequently the critical 
density for the transition, so we fix = 1 once for all). 

So far we defined the network and the dynamics along a link (street); we must further 
specify what happens when more links merge in a node (crossroad): several decision rules 
can be implemented; in this preliminary work we impose a simple random walk at the 
nodes: if a car is at the end of a link and has n possible directions to take it chooses one 
of them with probability n~^. 

Of course the total number of machines is conserved along the dynamics and we impose 
Kirchoff rules [2] for the flow at the nodes. 



3 The adapted Lighthill-Whitham theory 

As we want to move from a microscopic prescription toward a macroscopic description, 
there should be a mesoscopic lengthscale at which well known models, the most famous 
being the Lighthill-Whitham model (LW) [13], should be recovered. With mesoscopic we 
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mean we are dealing with an ensemble of cars such that their concentration as a func- 
tion of the space (labeled by x) and time (labeled by t) is meaningful whilst the average 
overall the concentrations can still be thought of in terms of a probabilistic description 



From a practical viewpoint let us now concentrate just on a big street with an amount 
M of cars inside and analyze the flow at this level. 

Assuming that the cars move from smaller to larger values of x, we define the concen- 
tration C of an element of the traffic on the street as the amount of cars moving within 
the space delimited by two generic points C D x E [xa,Xb], Xa < Xb- It follows that 
C{x,t) = Yld ~ ^a) — — Xb)) and let us also write the velocity of the 

generic z*^ particle as Vi{t) = 1 — 6{{xi{t) + d) — Xi+i{t)), such that if the car is moving, 
its velocity is one, otherwise is zero. 

The traffic flow is deflned as J{x,t) = Vg{x,t)C{x,t) where Vg{x,t) is the group velocity 
and C{x,t) the concentration. Assuming conservation on the total amount of cars, the 
following continuity equation holds 



'-^^'-^-tplA^^')-tpLM (3.1) 

« j 

where p\n{x, t) = 9{xa — Xi{t)) — 9{xa — Xb — Xi{t)) and pluti^, t) = 1 — p\n{x, t) take into 
account car entries and exits on the street. 

The group velocity relates to the local velocity of the particles via [2] = Yld '^'i + 
jjC^c ^ Vi < 1 which in our case reads off as 

N 



^. = (^'.)-(^E^.> = l-<^' (3-2) 



where the brackets average out the Dirac deltas on some measure (note that Vg < (vi)). 
The scenario is as follows: as far as the system is completely un-frustrated (low density 
liquid) the group velocity corresponds to the single velocity; increasing density of cars 
(or decreasing the connectivity of the network) frustration arises and lowers the group 
velocity up to a threshold where the jamming transition starts and Vg goes to zero such 
that the only solution to eq. (13.11) is C{x,t) = C e [0, 1]. 

For the sake of simplicity let us consider a very long piece of a street without entries or 
exits such that the LW model in our case obeys the following partial differential equation 

'-^.i.- ,f-£^ ^ 0, (3.3) 

whose solutions are Galilean characteristics with slope 1 — and can be expressed by 
generic functions f{x ± (1 — 4>)t) (cinematic waves [63])- If we now focus on two adjacent 
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V C 

Figure 3.1: Left: Collision time versus the fluidity of the first concentration Ci in our 
LW model, theor stands for analytical prediction, exp for the simulation results. Right: 
Fundamental diagram for our LW model. 

elements, Ci{x, t), C2{x, t) and suppose that Ci(x, t) is in front of C2(a;, t) in the direction 
of motion, we see that as far as their corresponding fluidities are zero they both follow 
straight lines (on the x, t plane) with the same slope. If we now suppose that Ci changes 
its status (for example by internal rearrangements, that we impose in the simulation by 
freezing randomly spheres belonging to the first package) and becomes frustrated (pi > 0, 
its group velocity decreases and at the time AX/{vgi{(f>) — Vg2) the two characteristics 
meet causing a discontinuity in the concentration functions, which acts as the starter of 
the transition, as (due to the structure of the solution) it propagates both forward and 
backward on the street. 



4 Macroscopic behavior 

Now we turn to the analysis on large scale networks. We consider three sizes of (squared) 
networks, which, at the highest connectivity, are completely filled by Nmax = 2 x 10'^, 

5 X 10'^ and 1 x 10^ spheres, which set p = 1 in the different cities. 

The connectivity ranges from to 1 and is changed from the percolation threshold (which 
is defined by a = 0) to the fully two dimensional network (which is defined by a = 1) in 
eight steps and the density (again defined in [0, 1]) is investigated at five different values 
too. 

The analysis works as follows: for every investigated size, degree of connectivity and den- 
sity, we average over twenty different runs; for every run we start off spheres randomly 
over the network (avoiding pathological overlaps) and study the mean of the fluidity, its 
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Figure 4.1: Fluidity fluctuations in the steady state for the largest size city at a connectiv- 
ity a = O.6.- We show the distributions of 3 different densities (N = 1 ■ 10^, 2 • 10^, 3 • 10^^ 
on a lin-lin scale (left), the spread of their variances, centering the distributions on the 
smallest average for easier visualization (center) and a Gaussian fit (x^ ~ 0.93) on a 
log-scale obtained by sampling over 2000 cars. 

variance and its distribution by starting to collect data when the mean does not vary 
sensibly (less than 5% for an amount of time O(IO^) collisions), that we claim to be a 
stationary state. 

At first, even though we do not provide a plot (as it would have just three points), 
we stress that we obtained good agreement among the three different city sizes investi- 
gating the monotonicity of the convergence of the averages (of fluidity and its variances), 
which is the basis of the thermodynamic limit (fundamental property for a well behaved 
model [5^). 

For all the not-jammed stationary states (at fixed a,p) we analyzed the distribution of 
the fiuidity sampling over the whole simulations performed: they turn out to be close 
to Gaussians distributions (which we use as a test fit, see figure (14. ip ). the variance be- 
ing almost independent by a and increasing with the density of the spheres up to the 
transition point (at given a), immediately later they are delta-like on the averaged of 
the fiuidity (0 = 1). The average fiuidity versus the density (that we plot directly by 
using the amount of cars) is shown in figure fl4.2p : For low density regimes, fiow behavior 
appears independent by the connectivity of the network (and the scaling is always the 
same (/> oc p) , while, for higher values of the density a second regime is approached which 
is highly sensible by the connectivity, and in which, continuity of the order parameter 
with respect to the density is still observed. 

Finally, for highest level of density a discontinuous jump to the jammed phase is ob- 
served, for all the values of the connectivity (see figure (14. 3 p where we report results for 
the medium city size). 

In figure fl4.3p furthermore we show the fiuctuations of the order parameter for the small- 
est city: it is worth noting that a phase transition (marked by a sharp peak inversely 
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Figure 4.2: Left: Behavior of the fluidity versus the total number of cars at six different 
values of the connectivity. In the low density regime the dependence of versus N 
is linear and universal. Right: Example of fluidity versus time on a single network 
(time is measured by the number of collisions (modulus 10^ j showed from the larger city: 
there exist only two "macroscopic" long term behavior: a jammed phase, that for the 
investigated values appears for N > 7 x 10^ and a liquid phase, which however shows a 
continuous variation in the fluidity versus the density of cars (for N < 7 x 10^). 




Figure 4.3: Left: Medium city. Behavior of the fluidity versus the connectivity for five 
different densities, shown in terms of the total amount of cars N . Right: Smallest city, 
self- averaging C of the order parameter: fluctuations versus the density ( expressed in 
terms of the total amount of cars), for six different connectivity values. The onset of the 
jam transition is shown by the peak, which, increases proportionally to . 
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Figure 4.4: Phase diagram: The line (p = 1 is the transition line. At right there is 
the jammed phase at left the fluid phase. Further, the latter is shown split in six zones 
accordingly with different values of the fluidity. 

proportional to the connectivity) seems to appear at a critical value of density. Analyzing 
again the medium size city, we show in figure (14.31) the behavior of the order parame- 
ter versus the connectivity (for several values of density) and there is no presence of a 
continuous behavior: at a critical value (depending on the density) a jump to a jammed 
phase is observed. 



5 Conclusions 



In this paper we implemented a numerical algorithm which mimics the fiow of cars in 
urban cities: cars are described as kinetically constrainedhsnd spheres and urban topology 
is chosen as a planar Erdos Renyi graph. Kinetically constrained because the one way 
roads break detailed balance and the collisions among two cars do not respect the third 
law of thermodynamics. Even though mathematically hard to be analyzed, this model 
can still be investigated by numerical simulations. Hard spheres because, as in principle 
we do not know how many timescales are involved in the genesis of the congested phase, 
and being interest in the long time behavior, we have chosen one of the fastest possible 
algorithms for the dynamics: the event driven motion, which requires hard spheres. We 
investigated the response of a dynamical order parameter, the fiuidity, defined as the 
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ratio among the moving cars on the whole ensemble, by tuning two control parameters: 
the density of the cars and the connectivity of the network. 

From this numerical investigation we found a continuous transition from a congested 
phase to a fluid phase by varying the density of the cars (at fixed connectivity) such that 
the fluidity lowers smoothly from 1 to smaller values (up to zero where there are no longer 
cars) and a discontinuous jump of this order parameter when varying the connectivity of 
the network at fixed amount of cars. 

Furthermore the timescales involved seem to be several (the longest of which seems to 
diverge at the transition to a jammed phase) rising the question on what kind of traffic 
optimizer should be developed in order to minimize traffic. 

On these first heuristic considerations we believe that interacting local optimizers (as 
a grid of interacting cross-lights able to detect flow [571 ES]) would work better than a 
global ground state searcher and are more stable with respect to perturbations as new 
added (or removed) streets. 

Future works concerning the kind of transition will be due to investigate the relaxation 
to equilibrium after the stimuli by introducing a car (or a few) or by introducing a new 
link, so to check the presence of aging in the network. Furthermore traffic optimization 
by properly interacting traffic lights will be considered as well. 
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